A space-time integral estimate for a large data semi-linear wave 
equation on the Schwarzschild manifold. 



P. Blue A. Soffer 

February 2, 2008 



Abstract 

Weconsider the wave equation (-d' 2 + d 2 p - V- V L (-A S 2))u = fF'(\u\ 2 )u with (t,p,6,<f>) inExlxS 2 . 
The wave equation on a spherically symmetric manifold with a single closed geodesic surface or on the 
exterior of the Schwarzschild manifold can be reduced to this form. Using a smoothed Morawetz estimate 
which does not require a spherical harmonic decomposition, we show that there is decay in L 2 oc for initial 
data in the energy class, even if the initial data is large. This requires certain conditions on the potentials 
V, Vl and /. We show that a key condition on the weight in the smoothed Morawetz estimate can be 
reduced to an ODE condition, which is verified numerically. 



We consider the following defocusing wave equation: 

(-d 2 + d 2 p -V- V L (-A s2 ))u =f(p)F\\u\ 2 )u (1) 

■u(O) =u 
d t u(0) =ui 

for a function u : R x M — > R with {t,p,6,4>) £ R x M. = R x R x S 2 . Our goal is to show that, under 
conditions fTTTTTI on the potentials (see below), 







if the initial data (ito,Mi) has finite (but possibly large) energy. In fact, we show the stronger result that 
So I M (l + P 2 y 2 H 2 dtdpd 2 Lu < CE[u Q , Ul }. 

Equation ([1]) describes several interesting geometric wave equations. On a spherically-symmetric, 3- 
dimensional, warped-product of R and S 2 , Riemannian manifold, the metric is ds 2 = dp 2 + r(p) 2 duj 2 , and if 
r has a unique minimum, it corresponds to a closed geodesic surface. This is the first problem we consider, 
and, by translation, we may assume that the closed geodesic surface occurs at p = 0. The semi-linear 
wave equation (— d 2 + A g )u = |{t| p_1 -u can be reduced to |T]) by setting u = ru, in which case V — r"/r, 
Vl = r~ 2 , and / = r 1_p = vjf 1 ^ 2 . The conditions [TlfTTI on the potentials are, therefore, conditions on r. 
These conditions are not vacuous, since, for example r = 1 + p 2 generates a wave equation which satisfies 
these conditions. Similarly the wave equation on the exterior region of the 3 + 1-dimcnsional, Schwarzschild 
manifold also can be reduced to ([T]) if p = 3, and this is the second case we consider. Both cases are explained 
in more detail in [2]. 

The purpose of this note is to show that there is decay even for large data (defocusing) semi-linear 
problems. In [2], we considered the linear equation (([T} in which the right-hand side is absent) in both the 
Riemannian and Schwarzschild cases, and the non-linear problem with small initial data in the Riemannian 
case. In these cases, we were able to show that the space-time integral || "«|| i 4 (t/tdpd 2 uj) is controlled by weighted 
norms. This built on previous work [3], in which the L°° norm was shown decay if weighted H 3 norms 
were bounded (and the initial data was small in the non- linear case). Similar results, with better decay 
estimates at the event horizon (an important part of the Schwarzschild manifold, corresponding to p — > — oo) 
have also been proven [3]. 

We make the following assumptions on the potentials and non-linear terms: 
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1. [Radial potentials] V and Vl are smooth functions of p only. 

2. [Positive energy] V and Vl are positive. 

3. [Unique, unstable, trapped geodesic surface] Vl (which we refer to as the angular potential) has a 
unique critical point, which is a maximum, and occurs at p = 0. We refer to p = as the trapped ray, 
trapped geodesic surface, or the orbiting null-geodesic surface. 

4. [Derivative of V well-behaved] d p V is bounded and goes to zero at ±oo. 

■ ,E+1 

5. [Semi-linearity] There's ap€ (1,3] such that F(x) = p+1 ■ 

6. [Compatible semi-linear coefficient and angular potential] / = V L 2 . 

Throughout this note, we use C to denote an arbitrary constant which may change from line to line, (•, ■) as 
the L 2 inner product (with respect to dxd 2 u> unless otherwise stated), and || • || as the L 2 norm. 

1 Analysis using the method of multipliers 

It is well-known that there is a conserved energy, as in R™ +1 , 

E[u,d t u}=[ \3 t u\ 2 + \d p u\ 2 + V\u\ 2 + V L \V S 2u\ 2 + 2f{p)F{\u\ 2 )dpd 2 Lo. (2) 

J{t}xM 

Following our earlier work [2], we introduce a radial derivative operator, pointing away from the trapped 

ray 

l = gd p +g'/2 
in terms of a C 3 radial function g — g(p) satisfying 

7. [Increasing] g'(p) = d p g > 0, 

8. [Centered about the trapped ray] g(0) = 0, 

9. [Energy bounded] g is bounded, 

10. [Inverse polynomial behaviour] there's a power a € [—2, —1), and for i 6 {1,2, 3}, there are constants 
Cj and Ci such that (for |p| sufficiently large) Cip a+1 ~ l < g^ < Cip a+1 ~ l , and 

11. [Positive spectral condition] There's an t\ > and a non- negative function xi which is positive in a 
non-empty open set, such that Wtp £ C°°({£} x M) : 

/ 2. 9 '|a p ^| 2 + (-g"'/2 - gV'm 2 dpd 2 Lo > ex [ g'\d p ^\ 2 + X M?dpd 2 Lo. 

J{t}xM J{t}xM 

This is equivalent to saying that the (self-adjoint) operator A = —(2 — ei)d p g'd p — g"'/2 — gV — eiXi 
has non-negative spectrum. 

Because condition ITTI involves V, the existence of g satisfying conditions [7lfTT1 is a condition on V. In the 
Riemannian case, with r = 1 + p 2 , since V — 2/r 2 has a unique maximum at p = 0, the arguments in [3] or 
[2] shows the existence of a g satisfying conditions [7]fTT1 (technically, the functions in [3] and [2] have g'" as 
a distribution, not a continuous function, but this is sufficient for the analysis). 
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Applying the method of multipliers to ju, we have 



((u,7«) + (•yu,u))\f 1 = 



2g'\d p u\ 2 + (-$"72 - gV')\u\ 2 + (~gV[)\\7 S 2 U \ 2 dfxPwdt 




+ 



/ «7u, /F'(| M | 2 )u) + </F'(| U | 2 ) U , t«)) eft . 



Term /V 



We now show that the integrand on the right is positive. By condition 1111 terms I and II are positive 
when taken together, and, in fact, dominate (1 + p 2 )~ 2 \u\ 2 through the Hardy estimate in lemma [T] The 
contribution from term III is positive because by conditions [3] and [5J g changes from negative to positive at 
the same point —V L does. At this point, we analyse the structure of the non-linear contribution, term IV in 
greater detail than in [SJ. 



Since Vl is positive and d p g > 0, the first term in brackets is positive. Since g changes sign at the maximum 
of Vl, —gdpVh > is also positive. Thus the integrand is positive. 

We now turn to the left-hand side. Since g is bounded, ||gd p u|| 2 < CE. Since g' < (1 + |p|) _1 , by the 
Hardy estimate, ||(9 p g)u|| 2 < CE. Thus, the left hand-side is bounded by the energy, and 



Taking t\ = and t% — > oo gives the desired result . 

We now prove a smooth Hardy estimate. Essentially, since Ai = K x S 2 is a 3-dimcnsional manifold, 
one expects the Hardy estimate from R 3 , ||u/|x||| < ||Vu|| to hold; however, because there is no origin for 
K x S* 2 , a little extra control is required. This result is taken from |3]. 

Lemma 1 (Smooth Hardy). For any a > and non-negative function \ which is positive in a non-empty, 
open set, there 's a constant C , such that 



( 7U JF'(| M | 2 ) M ) + (/F'(| U | 2 )u, 7M ) 





(3) 
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Proof. For p\ > 0, 

(1 + Pi) a+1 1 Wl Jo P (l + p) Q + 1 



Pl 2 U ^ p u M 2 
(1 + P) Q+1 V ; (l + p)«+ 2 

a +i r n 2 , . 2 r_MLw r r n 2 , 



2 Jo (1 + P)"+ 2 a+lj (1 + P ) a y 'J (l + P ) a + 



(l+p)«+ 2 r -(a + l) 2 J Q (l + p) a a + Wl 

We take the limit p\ — > 00. Since (for any exponent /? > 0) (1 + p)' 3 is equivalent to (1 + p 2 )' 3 / 2 on [0, 00), 
the powers of (1 + p) can be replaced by (1 + p 2 ) 1 ^ 2 . By symmetry, the same result holds on (—00, 0]. Since 
(1 + p 2 ) - ' 3 is uniformly equivalent to (1 + (p — po) 2 ) ^ for po in a finite interval, the |u(0)| 2 term can be 
replaced by |it(po)| 2 in any fixed interval. By integrating the estimate over po with po in a bounded open 
set in which \ is positive, the desired result holds. □ 

2 Using numerical ODE solutions to verify the spectral condition 

In this section, we (i) show that the spectral condition can be reduced to showing that there are no zeroes 
for the solution to an associated ODE and (ii) numerically verify this condition for some choice of function 
g on the Schwarzschild manifold. 

2.1 Reducing the spectral condition to an ODE condition 

Lemma 2. The spectral condition, follows from the following conditions: 
12. [ODE condition] There is a smooth, non-negative solution ipo to 

(-(2 - e 1 )d p g'd p - g'"/2 - gV - e lX i)V>o = 0. (4) 

Proof. We introduce several new variables to make this problem simpler. First, we change variables to 
x = x(p) defined implicitly by 

^ =<?', x(0) =0. 

ax 

Since g' is positive, x is well defined. By condition 1 10| x grows polynomially in p. Let 

A = - (2 - e 1 )d p g'd p - g"'/2 - gV - e lX i 
B = -(2-e 1 )d 2 x + W 

W =g'{ P {x)) (- g/,,( ^ (x)) - g(p{x))V\p{x)) - e lX Mx))) ■ 

Note that, by conditions HI [H andPTOl W — > as x — > ±00. If tjj e C°°(p,cj), then for <f> <E C°°{x,uj) given 
by <j>(x,w) = ip(p(x),uj), 
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j ipAipdpd 2 uj = J <pB(f>d. 



xd 2 uj. 
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Thus, the spectrum of A and B are the same, and if Aip = with -ip nowhere zero, then Bcf> = with <j) 
nowhere zero. Let </>o be a non-negative solution to B(j)Q — 0. 

If 4> € Co°' then, since 4>o is positive, we can write it as = 4>qu with u € . Consider the expectation 
value of B with respect to <p. In this paragraph, we use c = —(2 — ei) < 0. 

(</>, B<t>) = J (t) u{cd 2 x + W)<j} udxd 2 u; 

= J (j)QU ((cd 2 cj)o)u + 2c(d x (j)o)(d x u) + cj)o(cd 2 u) + W(j)ou) dxd 2 uj 

= J <t) u((B(j) )u + 2c{d x (j)Q){d x u) + (j) {cdlii))dxd 2 uj 

J (f>oU(f>o(cdlu)dxd 2 uj — J 4>q (d x (u(cd x ii)) — (d x ucd x u)) dxd 2 u> 

= J -2(t>o(d x (f>o)u(cdxu) - (j)1{d x ucd x u)dxd 2 oj 




This proves that B and, hence, A have non-negative spectrum. 

□ 

2.2 Numerical verification of the ODE condition on the Schwarzschild manifold 

To verify condition Q21 one can solve the ODE numerically. There are two potential problems with this 
method: first, the accuracy of the numerical solution, and second, solutions can only be found on a finite 
range. Because of the power of computers, we will expect that numerical solutions are accurate and ignore 
the first problem. We demonstrate how to control the errors caused by the second problem by considering 
the wave equation on the Schwarzschild manifold. 

On the Schwarzschild manifold, there is another radial coordinate r S (2M, oo) defined implicitly, and 
used to define the potentials, 

dr , 2M, . . 2M , 2AT 1 , 2M 

— =1 ), rO =3Af, V=— 1 ), V L =-(1 . 

dp r r° r r z r 

From the definition of r, it follows that 1 — 2M/r decays exponentially in p as p — > — oo and r/p — > 1 as 
p — > oo. Let 

The function g' decays like b~ 1 p~ 2 , g'" decays like 66~ 1 p~ 4 , and gV decays exponentially as p — > — oo and 
like Cb~ 1 / 2 p~ 4 as p — > oo. 

We will solve the differential equation Arp = numerically in the region (—po,Pa) and estimate the 
solution in the asymptotic region \p\ > p - 

We proceed with the asymptotic analysis first, with the goal of finding conditions on the numerical 
solutions to guarantee that the solution will remain positive in the asymptotic region. Since — gV is positive 
for large p, we will ignore it in the asymptotic region, leaving this part available to be used as xi- 

To investigate asymptotic behaviour, we use the Schrodinger operator B in terms of the variable x. We 
analyse the situation for p 3> 1, but the case p <C — 1 is similar. In this paragraph (and only this paragraph), 
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we use a ~ b to mean there is an e such that (1 — e)a < b < (1 + e)a, with e < 2/(1 + bp 2 ,). 



dx bp 2 



dx 

bp 3 ~3x 
6 
"bp 



m 6 
9 ~T7i 



-99 /2~-^62 
3 

~ ~ 9x2 

Thus, taking — (1 + e)/3x 2 as a lower bound for —g'g"'/2, it is sufficient to show that 

B 1 =-(2-e)%-(l + e)± 

has a positive solution to Bip = 0. Using e' = 3e < 1/100, we can replace this condition by positivity of the 
corresponding function for 

2 1 + 



B = -^ 3x 2 ■ 
Taking the ansatz <p — x a , we can find solutions 



-2 ± - 8 + 8e73 1,1^.,, 

a = = - ± — — + (J(e ) corrections. 

-4 2 2V3 

A condition must be found to fit the asymptotics to the numerical solution. In the region x — > oo, wc 
require the positivity of the coefficient, C\, on the more rapidly growing monomial, so that the solution will 
remain positive. Using 

cp(x) =dx ai + C 2 x a2 

<t>'{x) ^daxx^- 1 + C 2 a 2 x a2 - 1 

n <f>(x) - C\x^ 
= 



</>'{x) =Ciaix ai " 1 + a 2 (4>{x)x- x - Cix" 1 " 1 ) 
_ (/>'(x) - a 2 (f>{x)x- 1 ^_ ai+1 



ati — a 2 

Thus the condition we require is that, at the point x where we match the numerics to the asymptotics, 

d x (j> >a 2 (f>(x)x^ 1 . 

Thus, we require(again with e < 2/(1 + bp 2 ,)) 

1 3 
g'd p tp(p) >Y~^ a2 ^^bp3 

1 3 

dp^(p) >Y—^a 2 iP{p)-. (5) 
A similar result is required as p — > —00, but with the signs reversed. 
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We now show numerically that there is a positive solution to 



-(2-e 1 )d p g , d p -g , "/2-gV 

on (-po, Po) with 

p =1000, ei=l/1000, 6=.l 

(we also take M = 1 as a normalisation). This solution we construct satisfies 

^(-po)>-2(i-^M-po)^, (6) 

which, given the small size of e = 2/(1 + bp 2 ) ~ 1/50,000 and the additional factor of 2 relative to |(SJ), is 
sufficient to guarantee the matching of the numerics to the asymptotic solutions. 

We do this by treating Aip — as an initial value problem for ip(p) posed at p = po with initial conditions 
tf;(p ) = 1 and ip'(po) — 2(1/2— l/(2y/3))ip(po)-, finding the solution is positive, and verifying the condition 
is satisfied at —po- We show four plots: the potential —g"'/2 — gV', the solution in the left region [0,po]j 
the solution in the middle region [—10, 15] in which the potentials are large and the solution oscillates, and 
the solution in the right region [— pq,0] in which the solution goes rapidly to oo as p — ► — oo. Note that, 
although the asymptotic solution is concave down (with exponent less than 1) in terms of the variable x, 
since x — bp 3 /3, the solution is concave up in terms of p (and this behaviour is already clear from the plots as 
p — > ±po)- From the numerics, we find that at — po = —1000, the solution has value ~ 150000 and derivative 
~ —370. From ([6]), the derivative must be less than ~ —170. Since this is satisfied, the solution continued 
to the left of the numerical approximation will always be positive. By our choice of initial conditions, to the 
right of the numerical approximation, the continuation of the solution will also be positive. Thus, we have 
verified the ODE condition, condition \12\ and hence the L 2 oc result in ([3|). 
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